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Abstract 

The phenomenological motivations, the expressions and the comparison with 
experiment of the parton distributions inspired by the quantum statistics are 
described. The Fermi-Dirac expressions for the quarks and their antiparticles 
automatically account for the correlation between the shape and the first 
moments of the valence partons, as well as the flavor and spin asymmetries 
of the sea. One is able to describe with a small number of parameters both 
unpolarized and polarized structure functions. 



Let us first recall some of the basic ingredients for building up the parton 
distribution functions (PDF) in the statistical approach, as oppose to the 
standard polynomial type parametrizations, based on Regge theory at low x 
and counting rules at large x. The fermion distributions are given by the sum 

1 Invited talk at XX International Workshop on '"Deep-Inelastic Scattering and related 
subjects'", University of Bonn, Germany, March 26-30 2012, presented by F. Buccella, to 
appear in the conference proceedings 



of two terms [I], the first one, a quasi Fermi-Dirac function and the second 
one, a flavor and helicity independent diffractive contribution equal for light 
quarks. So we have, at the input energy scale Ql = 4GeV 2 , 

h, rfl\ = AX ^ xb ~ Axh 

xq (x, Q ) - ^ x _ X hjiz\ + 1 + exp(x/x) + 1 ' { } 

xq h (x, Ql) = q \ + — . (2) 

qy Wo) exp[(x + X^ g h )/x] + 1 exp(x/x) + 1 y } 

Notice the change of sign of the potentials and helicity for the antiquarks. 
The parameter x plays the role of a universal temperature and X^ g are the 
two thermodynamical potentials of the quark q, with helicity h — ±. It 
is important to remark that the diffractive contribution occurs only in the 
unpolarized distributions q(x) = q+(x) and it is absent in the valence 

q v (x) — q(x) — q(x) and in the helicity distributions Aq_(x) = q + (x) — q-(x) 
(similarly for antiquarks). The eight free parameters^! in Eqs. ( 11121) were 
determined at the input scale from the comparison with a selected set of 
very precise unpolarized and polarized Deep Inelastic Scattering (DIS) data 
P]. They have the following values 

x = 0.09907, b = 0.40962, b = -0.25347, A = 0.08318, (3) 

X+ = 0.46128, Xq u = 0.29766, X~ d = 0.30174, X+ = 0.22775 . (4) 
For the gluons we consider the black-body inspired expression 

xG(x,Ql) = j^jz-s r, (5) 

exp(x/x) — 1 

a quasi Bose-Einstein function, with be = 0.90, the only free parameter 
since Aq = 20.53 is determined by the momentum sum rule. We also assume 
that, at the input energy scale, the polarized gluon distribution vanishes, 
so xAG(x, Ql) = 0. For the strange quark distributions, the simple choice 



2 A = 1.74938 and A = 1.90801 are fixed by the following normalization conditions 
u — u = 2, d — d = X. 

3 In Ref. pQ we were assuming that, for very small x, xG(x, Qq) has the same behavior 
as xq(x, Qq), so we took be = 1 + b. However this choice leads to a too much rapid rise 
of the gluon distribution, compared to its recent determination from HERA data, which 
requires bij — 0.90. 
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made in Ref. pQ was greatly improved in Ref. [2J. More recently, new tests 
against experimental (unpolarized and polarized) data turned out to be very 
satisfactory, in particular in hadronic collisions, as reported in Refs. [31 H]. 

For illustration, we will just give one recent result, directly related to the 
determination of the quark distributions from unpolarized DIS. We display 
on Fig. 1 (Left), the resulting unpolarized statistical PDF versus x at Q 2 =10 
GeV 2 , where xu v is the w-quark valence, xd v the ci-quark valence, with their 
characteristic maximum around x = 0.3, xG the gluon and xS stands for 
twice the total antiquark contributions, i.e. xS(x) = 2x(u(x) + d(x) +s(x)) + 
c(x)). Note that xG and xS are downscaled by a factor 0.05. They can be 
compared with the parton distributions as determined by the HERAPDF1.5 
QCD NLO fit, shown also in Fig. 1 (Right), and there is a good agreement. 
The results are based on recent ep collider data from HERA, combined with 
previously published data and the accuracy is typically in the range of 1.3 - 
2%. 



HI and ZEUS HERA I+II Combined PDF Fit 
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Figure 1: Left : BBS predictions for various statistical unpolarized parton 
distributions versus x at Q 2 = lOGeV 2 . Right : Parton distributions at Q 2 = 
lOGeV 2 , as determined by the HERAPDF fit, with different uncertainties 
(Taken from Ref. [5]). 

Another interesting point concerns the behavior of the ratio d(x)/u(x), 
which depends on the mathematical properties of the ratio of two Fermi-Dirac 
factors, outside the region dominated by the diffractive contribution. So for 
x > 0.1, this ratio is expected to decrease faster for X^ d — x < x < Xq u + x 
and then above, for x > 0.6 it flattens out. This change of slope is clearly 
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visible in Fig. [2] (Left), with a very little Q 2 dependence. Note that our pre- 




Figure 2: Left : The ratio d(x)/u(x) as function of x for Q 2 = 4GeV 2 
(solid line) and Q 2 = lOOGeV 2 (dashed-dotted line). Right : Comparison 
of the data on d/u versus x, from E866/NuSea at Q 2 = 54GeV 2 [7J, with 
the prediction of the statistical model (solid curve) and the set 1 of the 
parametrization proposed in Ref. [H] (dashed curve). 

diction for the large x behavior, differs from most of the current literature, 
namely d(x)/u(x) — > for x — > 1, but we find d(x)/u(x) — > 0.16 near the 
value 1/5, a prediction originally formulated in Ref. [6]. This is a very chal- 
lenging question, since the very high-x region remains poorly known. 
To continue our tests of the unpolarized parton distributions, we must come 
back to the important question of the flavor asymmetry of the light anti- 
quarks. Our determination of u(x, Q 2 ) and d(x, Q 2 ) is perfectly consistent 
with the violation of the Gottfried sum rule, for which we found the value 
I G = 0.2493 for Q 2 = 4GeV 2 . Nevertheless there remains an open problem 
with the x distribution of the ratio d/u for x > 0.2. According to the Pauli 
principle, this ratio is expected to remain above 1 for any value of x. However, 
the E866/NuSea Collaboration [7] has released the final results corresponding 
to the analysis of their full data set of Drell-Yan yields from an 800 GeV/c 
proton beam on hydrogen and deuterium targets and, for Q 2 = 54GeV 2 , they 
obtain the ratio d/u shown in Fig. |5] (Right). Although the errors are rather 
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large in the high-x region, the statistical approach disagrees with the trend of 
the data. Clearly by increasing the number of free parameters, it is possible 
to build up a scenario which leads to the drop off of this ratio for x > 0.2. 
For example this was achieved in Ref. as shown by the dashed curve in 
Fig. [2] (Right). There is no such freedom in the statistical approach, since 
quark and antiquark distributions are strongly related. On the experimental 
side, there are now new opportunities for extending the d/u measurement to 
larger x up to x = 0.7, with the running E906 experiment at the 120 GeV 
Main Injector at Fermilab [10] and a proposed experiment at the new 30-50 
GeV proton accelerator at J- PARC [TTj . 
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Figure 3: Left : Quark and antiquark helicity distributions as a function of 
x for Q 2 = 3GeV 2 . Data from COMPASS |12j . The curves are predictions 
from the statistical approach. Right : Ratios (Aw + Art) /(u + u) and (Ad + 
Ad)/(d + d) as a function of x. Data from Hermes for Q 2 = 2.5GeV 2 [13] 
and a JLab Hall A experiment [H]. The curves are predictions from the 
statistical approach. 

Analogous considerations can be made for the corresponding helicity dis- 
tributions, whose most recent determinations are shown in Fig. [3] (Left). By 
using a similar argument as above, the ratio Au(x)/u(x) is predicted to have 
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a rather fast increase in the x range [Xq u — x, Xq u + x) and a smoother be- 
haviour above, while Ad(x) / d(x) , which is negative, has a fast decrease in the 
x range (X^ — x, X^ d + x) and a smooth one above. This is exactly the trends 
displayed in Fig. |3] [Right) and our predictions are in perfect agreement with 
the accurate high-x data. We note the behavior near x — 1, another typical 
property of the statistical approach, is also at variance with predictions of the 
current literature. The fact that Au(x) is more concentrated in the higher x 
region than Ad(x), accounts for the change of sign of g™ (x), which becomes 
positive for x > 0.5, as first observed at Jefferson Lab [II] . 

Concerning the light antiquark helicity distributions, the statistical ap- 
proach imposes a strong relationship to the corresponding quark helicity 
distributions. In particular, it predicts Au(x) > and Ad(x) < 0, with 
almost the same magnitude, in contrast with the simplifying assumption 
Au(x) = Ad(x), often adopted in the literature. The COMPASS data [15] 
give Au(x) + Ad(x) ~ 0, which implies either small or opposite values for 
Au(x) and Ad(x). Indeed Au(x) > and Ad(x) < 0, predicted by the 
statistical approach [1] (see Fig. [3] (Left), lead to a non negligible positive 
contribution of the sea to the Bjorken sum rule, an interesting consequence. 
For lack of space we only mention the extension to the transverse momentum 
dependence (TMD), an important aspect of the statistical PDF and we refer 
the reader to Ref. [16] . 

A new set of PDF was constructed in the framework of a statistical ap- 
proach of the nucleon. All unpolarized and polarized distributions depend 
upon nine free parameters for light quarks and gluon, with some physical 
meaning. New tests against experimental (unpolarized and polarized) data 
on DIS, Semi-inclusive DIS and also hadronic processes, are very satisfac- 
tory. It has a good predictive power, but some special features remain to be 
verified, specially in the high- a; region, a serious challenge for the future. 
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